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SUMMARY 


An analysis is presented to calculate damping in annular seals for an incompressi- 
ble fluid. Results show that damping in tapered seals (optimized for stiffness) is con- 
siderably less than that in straight seals for the same minimum clearance. Damping 
in rotating seals can promote fractional frequency whirl. 

Results were obtained with and without the fluid temporal acceleration term 3u/3t. 
Neglecting fluid acceleration led to errors in calculated damping of up to 16 percent. 
However, this is considered acceptable for a large class of practical application. 


INTRODUCTION 

It is now well recognized that annular pressure seals can have a considerable in- 
fluence on the dynamic behavior of rotors. This is not surprising when one considers 
that such a seal has the appearance of a journal bearing, although with larger clear- 
ances than is usual in bearing practice. Most of the force generated in a seal is due to 
the high velocity throughflow of sealed fluid. In a series of papers (refs. 1 to 3) Black 
and his coworkers calculated stiffness and damping in annular seals having constant 
clearance in the axial direction and sealing an incompressible fluid. Allaire, et al. 
(ref. 4) extended the analysis to seals with large eccentricities; however, they ne- 
glected the effect of fluid temporal acceleration without fully understanding the errors 
produced. In reference 5, Fleming showed that stiffness could be considerably in- 
creased by configuring the seal so that the clearance is greater at the inlet than at the 
exit. The higher stiffness can be beneficial in stabilizing rotors by shifting critical 
speeds; floating- ring seals can also benefit from higher stiffness with resulting longer 
life. 

It is recognized that damping, as well as stiffness, influences seal behavior. Thus 
the purpose of this report is to calculate damping for the tapered seal configurations 
which were analyzed for stiffness in reference 5 . A second purpose is to fully deter- 
mine the errors produced by neglecting fluid acceleration. If reasonable accuracy can 
be obtained with the simpler procedure, analyses of more complex situations (e.g. , 
compressible fluids) will be greatly eased. 


SYMBOLS 


A integration constant 

B damping, N- sec/m 

B dimensionless damping, BC^jl^D-^p^p 

C seal radial clearance (concentric), m 

D seal diameter, m 

E inertia coefficient, kg 

E dimensionless inertia, ECg/pDL 3 

e seal displacement, m 

F seal force, N 

f integral (eq. (26)) 

G H - 1 

H clearance ratio, C^/Cg 

h local film thickness, m 

J defined by eq. (18) 

K stiffness, N/m 

K dimensionless stiffness, KCg/p^LD 

L seal length, m 

M seal leakage flow, kg/sec 

M dimensionless leakage, M^iDCg-^J^p 

o 

p pressure, N/m 

R seal radius, m 

S defined by eq. (19) 

t time, sec 

u fluid velocity, m/sec 

V ah/at 

x, y coordinates normal to seal axis 
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z axial coordinate 

a taper angle 

A 9/0t 

e eccentricity ratio, e/C 2 

rj total entrance loss factor, 1 + | 

0 angular coordinate 

X friction factor 

£ entrance loss factor 

O 

p fluid density, kg/m 

a \L/C 2 

co angular speed, rad/sec 

Subscripts: 

p precessional 

r radial direction or rotational 

s steady concentric condition 

t perturbed condition or tangential direction 

0 upstream stagnation condition or zero order in time 

1 seal entrance or first order in time 

2 seal exit or second order in time 

ANALYSIS 

The configuration to be analyzed is that of an annular clearance seal whose clear- 
ance decreases in the flow direction (fig. 1). Around the circumference, the clearance 
is given by 

h(z, 0 , t) = C(z) + C 2 e(t) cos 0 

where C is the clearance for a concentric seal. This implies that the only motion con- 
sidered is in the plane 0=0. 

C(Z) = - Q!Z 
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The notation = C(0) and Cg = C(L) is used. When the seal is eccentric (e * 0) or 
has a velocity or acceleration in the radial direction, a radial force is generated. To 
determine this force, the pressure distribution within the seal will be calculated and 
integrated over the seal area. In common with references 1 to 5, the following assump- 
tions are employed: 

(1) The fluid is incompressible. 

(2) Rotational effects on the flow field are neglected; the flow is one-dimensional 
in the axial direction . 

In addition, as in references 1, 2, and 5: 

(3) Eccentricity is small compared with the concentric clearance, that is, e « 1. 

(4) Friction factor is constant eveiywhere within the seal. Reference 3 shows that 
this introduces some error for the Reynolds number range where the friction factor is 
Reynolds number dependent. In the case of a high Reynolds number, however, friction 
factor no longer varies with Reynolds number (ref. 6). 


Development of Equations 

As in reference 5, we start with the continuity and momentum equations. These 
expressions differ from those of references 1 and 4 because the clearance h varies 
along the length in a tapered seal. 



+ u— + 

dZ 


ah 

at 


= o 


_I9e = x^ + u-^ + -^ 

p 8z h dz at 


The boundary conditions are as follows: 


P = P 0 = n 

u 2 


at z = 0 


p = 0 at z = L 


( 1 ) 

( 2 ) 


( 3 ) 


Variables h, u, and p will now be expressed as the sum of a steady state concentric 
value and a perturbation which will be assumed small compared to the steady- state 
value. Thus, 
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I? 


h = C + h t ^ 

u = u s + "t ?• 

P = P S + Pt j 


( 4 ) 


The continuity and momentum equations may then be divided into separate expres- 
sions involving steady-state and perturbed variables. The steady-state expressions 
are 

C — — ^ + u — = 0 (5) 

dz s dz 


! dp 

p dz 


C 


- = X — + u 


du 
s 

dz 


with the following boundary conditions: 


-\ 


r?pu c 


p s = Pq- 


at z = 0 


P s = ° 


at z = L 


( 6 ) 


( 7 ) 


where 7} = 1 + £ is the total entrance loss factor. The perturbed variable expressions 
(in which products of two perturbed variables are neglected) are 


dz Sz T dz St 


p Sz 

The boundary conditions are 


_i!h= 2 ^_ x vi +u 

C c 2 s Sz * dz St 


Pj. = -T]pu g Uj. at z = 0 
Pj. = 0 at z = L 


( 8 ) 

( 9 ) 


( 10 ) 


Steady state equations (5) and (6) are easily solved: 
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( 11 ) 


u l c l 

U = -±-± 

s C 



u. 



The solution of equation (8) is 


u ,A_ h t u l C l. ) .V 
‘ C 0 2 a 


where A is an integration function and 


V 



( 12 ) 


(13) 


(14) 


(15) 


There are alternative methods for solving the perturbed momentum equation (9). One 
method is to retain all the terms and solve the equation as Black did in reference 1 . 
Another method is to neglect the fluid temporal acceleration term, 9u^./9t, as Allaire, 
et al. , did in reference 4. This would introduce greater approximation but the accu- 
racy may be acceptable given the much easier solution procedure . Both methods will 
be presented in order that a comparison of the results may be made. 


Solution Retaining Fluid Temporal Acceleration 3u/9t 


Equations (11) and (14) are substituted in equation (9) to obtain 


l^t 

p 9z 


3hj.ufC 


2 p 2 

1 C 1 


2 u iCi A 2X^0, V : VV 

(\ + a) + — ±-±- (\ + a) + + - + -¥- 


a C* 


a 


(16) 


where the dot denotes differentiation with respect to time; that is, A = 3A/9t and 
V = 3V/3t = 3 2 h^.^3t 2 . This expression may be immediately integrated and the boundary 
condition p^. = 0 at z = L applied to obtain 
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, 2 ^l c l v (l 

a 2 \ C 


xVA ln _C + V (z . L) 

C 2 / a C 2 a 


(17) 


The boundary condition at z = 0, equation (10), may be applied to find the integration 
function A. As a first step, all terms containing A are moved to the left side of the 
equation; thus, 


u i A [' 


rj + (ct + G)(H + 1)] + In H = u 2 h,. [t? + (cr + G)(H 2 + H + 1)] 

u, VHL 


where H is the film thickness ratio C^/Cg, G = H - 1, cr = AL/Cg, and a has been 
replaced by its equivalent in terms of film thickness ratio: a. = GCg/L. 

This equation cannot be solved in the conventional manner because both A and its 
time derivative appear. However, the equation may be rewritten using differential op- 
erators (ref. 7). 




J+ HLtaH A 


G 


') A = 


n U. HL ttt ^ 

ufs - — (rj + 2<j)A - — 


where A = d/dt and, to simplify the writing, the substitutions 

J = 17 + (<r + G)(H + 1) (18) 


and 


S = n + (a + G)(H 2 + H + 1) (19) 

have been made . This has the appearance of an algebraic equation and, indeed, it is 
common practice to treat A as an algebraic quantity . The following procedure to de- 
termine A was used by Black in reference 1 . A straightforward algebraic solution for 
A yields 
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( 20 ) 


u 2 S - 


A = 


u^HL(rj + 2<j) 


A-S^A 2 


J + HLtaH A 

-L n 


h 


For high axial flow rates, the first term of the denominator is much larger than the 
second; hence a binomial expansion of three terms will be a good approximation. 
Truncating at three terms is equivalent to neglecting the effect of time derivatives 
higher than the second on the seal pressures. With this substitution, equation (20) be- 
comes 


A 




UjHL 

G 


(tj + 2cr)A - 




HL In H ^ , H 2 L 2 In 2 H 
U 1 JG u 2 J 2 G 2 



( 21 ) 


This is next multiplied out, and coefficients of various powers of A collected, ignor- 
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ing coefficients of A and higher powers . Making the definition 

A = (A q + A x A + A 2 A 2 )hj. (22) 


we obtain 


A 


1 


A 2 ” 


HL 

UjJG 



= _HL/SJnH + T? + 
GJ \ J 

SHJs!H + HlnH ( 
J 2 G JG 



(23) 

(24) 

(25) 


With the axial pressure profile known, equation (17) may now be integrated to determine 

the force on an axial strip. (It is not necessary to consider the contribution of p , as 

s 

it is uniform around the circumference and the net force from it is zero . ) After some 
algebraic manipulation, 
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f = 


f L p t| - l2 , v 

J P 2C 2 c 2 g 2 gc 2 V G I 


y L 3 Au 1 L(a + G) 
’ 2GC^ C 2 

Similar to equation (22), define 

£ - (£ 0 + A + f 2 A 2 )ht 

Then, substituting equations (22) to (25) in equation (26) results in 

u?L 


£ ° ° cT (I7+G I ( H + 2 " j) 


( 26 ) 


(27) 


(28) 


f - U 1 L 2 ctH 

1 " gc 2 g 




S In H 


r\ — — | 

/H In H A 

) J 

V G /_ 


f 2 “ 


c 2 g 


JLfsj^ + ^ + 2 \/l L lnH. 1 \_l 
GJ \ J / \ G / 2 


(29) 


_ Hflg + G) |~SH In 2 H + H In H ^ + 2tr) - 1 
GJ 2 GJ 




(30) 


The total seal force is now found by integrating around the seal circumference 



cos 0 R d0 


The perturbed clearance hj. varies around the circumference according to 

ht = \ 0 cos 0 

where h^ 0 is the amplitude of the perturbation at 0 = 0. Thus 


(31) 


(32) 
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(33) 


F = J p(f Q + t 1 A + f 2 A 2 ^ cos 2 e R d e 
F = 7rpR(f 0 + f x A + f 2 A 2 )^ 


In terms of stiffness, damping, and inertia coefficients, 

F = (K+BA + E A 2 )!^ (34) 

It is convenient to present results in terms of the following dimensionless coefficients: 
Stiffness 


- KCr 
K = 


P 0 LD 


= — (a + G) ^— + H - — ^ 
J \2 J/ 


(35) 


Damping 


BC r 


B = 


L 2 Dy ^ GV2J L 


2gH /in H 1 

\ + — (CT + G) / S - — H + 1 ) + 2cr\ 


/h In H A 

_ G \ G 

/ j V j / 

J 

\ G /_ 


(36) 


Inertia 

_ EC f 
E * 


pDL' 


= /SJ^ + ^ + 2 \/HlnH_A 

2G GJ \ J / V G / 


_ i _ H(cr + G) j 
2 J 


SH In- H + HinH (T? + 2ff) 


GJ 


GJ 


-]] 


(37) 


It may be verified that these results reduce to those presented by Black and Jensen 
(ref. 2) as H tends to 1. 

Solution Neglecting Fluid Temporal Acceleration 9u/9t 

When the fluid temporal acceleration is neglected, the momentum equation (9) be- 
comes 
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i 


( 38 ) 


1 QPt u s u t U o h t 

p d z C r 2 s 9z “t dz 


Equations (11) and (14) for u g and u^., respectively, are substituted in this, the result 
integrated in z, and the boundary condition at z = L applied to obtain 

+ — C -' ^ (39) 

a \a /lc Cj 

The integration function A is readily found by applying the boundary condition (eq. (7)) 
at z = 0 resulting in 

A=± [ W - + ( 40) 

Substitution of this into equation (39) and integration to determine seal forces is 
straightforward. The result for seal stiffness is identical to that for the more com- 
plete solution, as expected. For damping, 


B = — in + 2( 

GV2J *- 

The inertia term does not appear. 


(- +1 ) 


V G / 

V G /J 


RESULTS 

The preceding analysis has been used to calculate values of seal damping and in- 
ertia terms for tapered seals. Straight seal (zero taper angle) results are also pre- 
sented for comparison; they were calculated from the expressions in references 1 
and 2. It also may be verified that the expressions for tapered seal results presented 
herein reduce to those for straight seals in reference 2 as the clearance ratio H tends 
toward 1 . Tapered and straight seals are compared using the same minimum clear- 
ance, since this is the significant dimension in practical seal design. 
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Values of dimensionless damping coefficient B appear in figure 2 for a range of 
dimensionless seal length cr and with entrance loss factor £ = 0. (For turbulent flow 
£ is at most 0.2; results are not much different than for 4=0. This was also shown 
by Black (ref. 2) for straight seals.) Results are shown for both straight seals and for 
seals whose taper is optimized to maximize the stiffness to leakage ratio (ref. 5). For 
each seal, results are shown for the full solution (retaining 8u/9t) and for the approxi- 
mate solution (neglecting 9u/8t). 

Damping values are always higher for a straight seal, as would be expected since 
the average clearance is lower. The difference is least near a length c r of about 1.5, 
where the straight seal damping is minimum. Tapered seal damping is nearly constant 
out to <7 ~ 3; damping then increases with a. The approximate solution errs in oppo- 
site directions for straight and tapered seals. For straight seals, the solution neglect- 
ing 0u/9t overestimates damping up to 16 percent. For the optimum tapered seals, 
the approximate solution underestimates damping a maximum of 10 percent. These 
errors are considerably greater than those implied by the authors of reference 4, but 
are still considered quite acceptable by the author. 

The inertia term E is shown in figure 3 for the full solution. (The approximate 
solution does not yield the inertia term at all. ) For the optimum tapered seal E is 
very nearly constant over the entire range of seal length shown. For the straight seal 
E is negative for very small seal length; it increases rapidly with length and for large 
lengths becomes nearly constant at a value somewhat greater than for a tapered seal. 
Negative results are not shown in reference 2, but this appears to be due simply to in- 
sufficient points having been plotted for small seal length. 


Effect of Seal Rotation 


The calculations so far have neglected any seal rotation. For a rotating seal, ac- 
cording to the coordinate transformation used by Black and Jenssen (ref. 2), the forces 
on the seal journal are given by 


F v 


K - A Ec4 

- Bu; 


— — 

X 


B 

Eco 


X 


p 

E 

0 


• • 

X 

X 

= - 

4 r 

2 r 



- 


r 



- 





_ F y_ 


- - Bco 

L 2 r 

K - i Eoj^ 


y 


-Eo> r 

B 


y 

— — 


0 

E 


•• 

y 


( 42 ) 
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In equation (42), u> r is the shaft rotational speed. Recall that rotational effects were 

neglected in determining the seal coefficients K, B, and E . Thus one would expect 

equation (42) to be valid only for small rotational speeds; that is, Rco r « u^. 

The seal journal is now considered to be precessing in a circular orbit about the 

seal center with orbital frequency u>p . For this condition we will calculate the radial 

restoring force and the tangential (whirl direction) force . Since the seal has circular 

symmetry, it is sufficient to examine the situation when the seal journal center lies 

** 2 

along the x-axis, that is, when x = e and y = 0. Then x = 0, y = ecOp, x = -eoT, 
y = 0. The radial and tangential forces are 


F 


r 



(43) 


These expressions may be rewritten as 



(44) 


Note that e is by definition in the radial direction. Two conclusions may be drawn 
from these expressions: 

(1) The inertia term is a perturbation on the seal radial restoring force. If only 
whirl ratios Wp/a> r from 0 to 1 are considered, this perturbation has a maximum 
value at 0 and 1 and becomes zero for w p/ w r = 1/2. 

(2) The tangential force promotes forward whirl when cu /w r < 1/2, inhibits whirl 
when a>p/o> r > 1/2, and is zero for half-frequency whirl (a>p/co r = 1/2). This is re- 
markably similar to the situation in a full circular journal bearing. The seal differs 
from a journal bearing in the radial force being a maximum for Wp/u r = l/2 (assuming 
E > 0, which is always the case for a tapered seal and is true for all but small seal 
lengths for a straight seal). We may also infer that damping in a rotating seal is not 
entirely beneficial. The damping term produces a tangential force which inhibits syn- 
chronous whirl but promotes whirl at half -frequency or less. 
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Numerical Example 


Seal damping and inertia have been presented in dimensionless form, as were stiff- 
ness and leakage flow in reference 5. Tapered and straight seal stiffness and leakage 
flow from reference 5, as well as optimum clearance ratios, appear in figures 4 to 6. 
Also, damping and inertia forces vary with rotor speed. Thus the relative size of the 
various terms can only be compared for specific examples. The example chosen here 
is the space shuttle main engine (SSME) high pressure fuel turbopump (HPFTP) inter- 
stage seals that were examined in reference 5. These seals are 76 millimeters in 
diameter, 43 millimeters long, and have a radial clearance of 0.13 millimeters. The 
fluid sealed is liquid hydrogen; the pressure difference across the seal is 140 bars at 
the rotor design speed of 37 000 rpm (ref. 4). As determined in reference 5, the 
dimensionless seal length <r is 2.4. Seal coefficients are presented in table I for 
straight seals and for optimum tapered seals . Tapered seal stiffness is more than 
twice that of the straight seal. The tapered seal damping and inertia terms are 20 to 
25 percent less than the corresponding straight seal quantities . Also shown in table I 
are the terms Bo> r and Eo>^ which appear in equation (44) . These terms have the 
dim ensions of stiffness and may be compared with the stiffness term. For the straight 
seal, damping in the form Bco r is much greater than the stiffness term, while for the 
tapered seal Bu> r is somewhat less than K. The inertia term Ew r is quite small 
for both straight and tapered seals . 

The coefficients from table I may be used with equation (44) to yield radial and 
tangential forces per unit radial displacement for various whirl ratios (ratios of pre- 
cessional to rotational speed) . These are shown in table II for whirl ratios of 0 (no 
precession), 1/2, and 1 (synchronous precession). For the case of no precession, 
radial forces are reduced slightly from the stiffness of table I due to the inertia term. 
The damping term transforms to a negative tangential force, which means that there is 
a whirl-promoting force on the seal journal. This tangential force is analogous to the 
tangential force in a journal bearing. For the straight seal, the tangential force is 
somewhat larger than the radial force, while for the tapered seal, tangential force is 
only one-third of the radial force . For a whirl ratio of one-half, the inertia effect dis- 
appears, and the radial force is equal to the nonrotating force of table I, while the tan- 
gential force is zero. For synchronous precession (o^/oj^, = 1), radial forces are 
identical to the zero-precession case, and tangential forces are identical in magnitude 
but of opposite sign to the zero-precession case. The tangential force is thus a strong 
inhibitor of synchronous precession. 
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SUMMARY OF RESULTS 


An analysis was performed to calculate the damping in tapered annular seals for 
an incompressible fluid. Results show that damping in optimized tapered seals is con- 
siderably less than that in straight seals for the same minimum clearance. It was also 
pointed out that damping in rotating seals is not entirely beneficial, as it can some- 
times promote whirl. 

Results were obtained with and without the inclusion of the fluid temporal acceler- 
ation term du/dt. Neglecting fluid acceleration led to errors in calculated damping of 
up to 16 percent . However, this is considered acceptable because of the much greater 
ease of obtaining a solution. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 16, 1979, 

505-04. 
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TABLE I. - STIFFNESS, DAMPING, AND INERTIA 
COEFFICIENTS FOR HPFTP STRAIGHT AND 
TAPERED SEALS; ENTRANCE LOSS £ = 0 


Straight seal | Tapered seal 


Dimensionless stiffness, K 0.11 

Stiffness, K, MN/m 38 


Dimensionless damping, B 
Damping, B, (kN)(sec)/m 
Ba; r> MN/m 

Dimensionless inertia, E 
Inertia, E, kg 
Ec o 2 v , MN/m 


TABLE II. - RADIAL AND TANGENTIAL FORCES 
PER UNIT RADIAL DISPLACEMENT FOR HPFTP 
STRAIGHT AND TAPERED SEALS 
[All forces in MN/m.] 


K r = 88 


Straight seal 

K r = 

36 

K t = 

-39 

K r = 

38 

Kt = 

0 

K r = 

36 

K t = 

39 


K t = -30 























p ~~ 



Figure 1. - Tapered annular seal. 



Seal length, o = XL/C 2 


Figure 2. - Damping in tapered and straight seals. En- 
trance loss E, = 0. 
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Figure 5. - Leakage flow through tapered and straight 
seals. Entrance loss E; = 0. 


E o> 

o 

=3 CZ 0 * 

.E 

■s.8 2 
o o 


CVI 

o 



.1 .2 .4 .6.8 1 2 4 6 8 10 20 

Seal length, a = M/C 2 


Figure 6. - Optimum clearance ratio for tapered seal. En- 
trance loss E, = 0. 


i 


1. Report No. 2. Government Accession No. 

NASA TP- 1646 

3. Recipient's Catalog No. 

4. Title and Subtitle 

DAMPING IN TAPERED ANNULAR SEALS FOR AN 
INCOMPRESSIBLE FLUID 

5. Report Date 
April 1980 

6. Performing Organization Code 

7. Author(s) 

David P. Fleming 

8. Performing Organization Report No. 

E-124 

10. Work Unit No. 

505-04 

9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Paper 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

An analysis is presented to calculate damping in annular seals for an incompressible fluid. Re- 
sults show that damping in tapered seals (optimized for stiffness) is considerably less than that 
in straight seals for the same minimum clearance. Damping in rotating seals can promote frac- 
tional frequency whirl. Results were obtained with and without the fluid temporal acceleration 
term 9u/3t. Neglecting fluid acceleration led to errors in calculated damping of up to 16 per- 
cent. However, this is considered acceptable for a large class of practical application. 


17. Key Words (Suggested by Author(s) ) 
Annular seals 


Ring seals 


STAR Category 37 

Seal damping 




19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

Unclassified 

Unclassified 

21 


18. Distribution Statement 

Unclassified - unlimited 


22. Price* 

AO 2 


* For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1980 


National Aeronautics and 
Space Administration 

Washington, D.C. 

20546 

Official Business . 

Penalty for Private Use, $300 


THIRD-CLASS BULK RATE 



3 1 ID, D, 021580 S00903D5 

DIPT OF THE AIR FORCE 
AF WEAPONS LABORATORY 


ATTN: TECHNICAL LIBRARY (SDL) 
KIRTLAND AFB NM 87117 








